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Abstract 

Starting from the axial heavy meson exchange currents, constructed earlier in conjunction with 
the Bethe-Salpeter equation, we first present the axial p-, u— and ai meson exchange Feynman 
amplitudes that satisfy the partial conservation of the axial current. Employing these amplitudes, 
we derive the corresponding weak axial heavy meson exchange currents in the leading order in 
the 1/M expansion (M is the nucleon mass), suitable for the nuclear physics calculations beyond 
the threshold energies and with wave functions obtained by solving the Schrodinger equation with 
one-boson exchange potentials. The constructed currents obey the nuclear form of the partial 
conservation of the axial current. We apply the space component of these currents in calculations 
of the cross sections for the disintegration of deuterons by low energy (anti)neutrinos. The deuteron 
and the final state nucleon-nucleon wave functions are derived (i) from a variant of the OBEPQB 
potential, and (ii) from the Nijmegen 93 and Nijmegen I nucleon-nucleon interaction. The extracted 
values of the constant L\ a, entering the axial exchange currents of the pionless effective field theory, 
are in a reasonable agreement with its value predicted by the dimensional analysis. 

PACS numbers: 11.40.Ha; 25.30.-c; 25.60.-t 
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I. INTRODUCTION 



Since quarks are confined, the quantum chromodynamics [l|, 0] is not directly suitable for 
the investigation of nuclear physics phenomena at low and intermediate energies. Instead of 
quarks and gluons one employs effective degrees of freedom (hadrons), and ideas and meth- 
ods based on the concept of the spontaneously broken chiral symmetry 0,0]. The hadronic 
degrees of freedom, relevant for describing a nucleus and its response to the external elec- 
troweak interactions, are nucleon, A(1236) isobar and low lying mesons at the hadron mass 
scale, like 7T-, p-, lo-,... mesons. The interaction of the vector mesons with the baryons and 
pions is fixed by the vector dominance model (VDM) J5 . This concept has been developed 

10, 

the Standard Nuclear Physics Approach (SNPA) [1^ 



successfully during the last 4 decades @, 0, H I, H Q, and is called by some authors as 
the Standard Nuclear Physics Approach (SNPA) [l 21 ] . 
Starting from the early 1970s, a particular effort was devoted to the study of mesonic degrees 



of freedom in nuclei by investigating meson exchange currents (MECs) effects [13J, [l4j, [15, 
16, 17, 18, 19, 20, 2lj]. One of the best proofs of presence of the pionic degrees of freedom 
in nuclei follows from the study of the transition + <-> 0~ in the A=16 nuclei [l6|, 21], 
induced by the time component of the axial current. Detailed studies have shown that the 
experimental data and the calculations can be reconciled only if the time component of the 
weak axial soft pion exchange current 



14 



23 



is taken into account. The hard pion corrections 

change the result by 10-15 %. 
The time component of the weak axial MECs plays an import ant jx)le also in interpreting 

and on the + 



=12 nuclei 24 



0" 



the data on the isovector + <-> 1 + transition in the A= 
transitions in medium-heavy and heavy nuclei. In these last transitions, the part of this 
component arising from the heavy meson exchanges, contributes sizeably (2^,0, 27]. The 



weak axial exchange charge densities, derived for the phenomena at the threshold, are given 
by pair terms that are related to nucleon-nucleon potentials [25j, or are obtained from the 



26 



chiral Lagrangian 
In the classification of Ref. 
of the pion range is of the 
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the leading term of the space component of the axial MECs 
order ~ 0(1/M 2 ), where M is the nucleon mass. Being of 
relativistic origin, it is model dependent. This component of the weak axial MECs plays an 
important role in such fundamental reactions as 

(i) the neutrino reactions in nuclei at low and intermediate energies, 

(ii) the weak transitions in light nuclei, in particular the tritium beta decay, ordinary muon 
capture in 2 H and 3 He *, and the solar fusion pp and p 3 He processes, so important for the 
determination of the flux of the solar neutrinos, 

(iii) and the parity violating (PV) electron scattering that aims at elucidating the strange 
quark contribution to the electromagnetic structure of the nucleon 

With the use of the chiral Lagrangians, the structure of the space component of the weak 
axial MECs of the pion range was studied in detail in the hard pion model J^J in Refs. (l8l . 

and applied to various reactions in the lightest nuclei (39|, Sj], |H, l43j . The largest 
effect arises from the A excitation current of the pion range, that is partially compensated 
by the analogous current of the p meson range. Besides, the pair current of the p meson 



* For the recent review on these two reacti 
t Several calculations of the PV inclusive 131 




3G 



37[ electron scattering off 



33l [3J, |35| | and exclusive 

deuterons have already been done considering different theoretical issues, however, the axial exchange 
currents have not yet been included in the calculations. 
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range was used during the last decade in Refs. 12. l44j|. However, the origin of its derivation 
remains unclear. 

The method used in Ref. 25j for the construction of the axial charge at the threshold was 



45j to the derivation of the space component of the weak axial potential 



applied in Ref. 
MECs. 

The weak axial MECs of the Nitauj system have recently been studied in Ref. 46j, employing 



the Lagrangian based on the linear o model [471 ] . The model suffers from some problems 



discussed in detail in Ref. 48 



The weak axial MECs of the pion range in the Bethe-Salpeter approach to the nuclear two- 
body problem has been studied in various chiral models in Ref. [49|. Another recent con- 
struction of the weak axial one-boson exchange currents for the Bethe-Salpeter equation has 
been done in Ref. [lo| making use of the chiral Lagrangians of the NAnpooai system 38|, 51 



The obtained current operators fulfill the Ward-Takahashi identities and the matrix ele- 
ment of the full current, sandwiched between the two-body solutions of the Bethe-Salpeter 
equation, satisfies the PCAC constraint. 

The chiral Lagrangians 38, 51] are constructed in such a way j3, S, 
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43, 151|, [52, 153|, |54j, 



that they reproduce results obtained from the current algebra and PCAC in the tree 
approximation. Besides possessing the chiral symmetry, our Lagrangians are characterized 
by the following properties: (i) They respect VDM, reproduce universality, KSFRI, KSFRII. 
(ii) They provide the correct anomalous magnetic moment of the a\ meson, (iii) They 
reproduce the current algebra prediction for the weak pion production amplitude. It was 
explicitly shown [5lJ for the Lagrangian, based on the hidden local symmetry approach 
that if the a\ meson field is eliminated, the resulting Lagrangian preserves the properties 
(i) and (iii). Moreover, if the p meson field is eliminated from this Lagrangian, the chiral 
Lagrangian of the nucleons and pions is recovered j57[. Hence our Lagrangians consistently 
combine the chiral approach with the VDM concept, and provide a reasonable approximation 
at the tree level to the hadron amplitudes up to the energy scale ~ 1 GeV (~ m p , m ai , M). 
Subsequently, these Lagrangians were applied to the construction of the weak axial MECs 
in the tree approximation: the generic relativistic Feynman tree-level amplitudes satisfy 
the PCAC constraint and the nuclear MECs derived from them are required to satisfy the 
nuclear PCAC constraint (see below). In this approach, the weak hadron form factors are 
naturally of the VDM form, but the strong nucleon form factors should be introduced by 
hand. 

In practical calculations, one makes the non-relativistic reduction of the currents by expand- 
ing in Q/Mh, where Q is the momentum of the external particles or the momentum transfer, 
and Mh is the heavy meson or nucleon mass. Besides the leading order terms, the leading 
relativistic corrections have been calculated in the electromagnetic sector for the MECs of 
the pion range [ 



Due to the presence of the small expansion parameter, the model 
currents are expected to be valid in the energy /momentum region up to ~ 0.4 GeV . Let 
us note that the cross sections for the backward deuteron electrodisintegration calculated in 



5l| describe well the data up to Q 2 « 1.2 (GeV/cf [60j, [61 



We shall call the above described approach, where the currents and potentials are constructed 
in the tree approximation, as Tree Approximation Approach (TAA). The advantage of this 
approach is a relative simplicity and transparency. On the other hand, the obtained results 
can be considered fully realistic. Moreover, the nuclear PCAC constraint connects a part of 
the axial nuclear MECs (potential MECs) with the nuclear potential derived within the same 
approach. Making use of this potential in the production of the nuclear wave functions, one 
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can do consistent calculations of the potential MECs effect. Let us note that the problem 
with the consistency of the calculations is by most authors overlooked [62| . 
In principle, the TAA is improved by effective field theories (EFTs). An EFT is based 
on the most general Lagrangian involving the relevant degrees of freedom, and respecting 
chiral symmetry [l], 0, [63[ . For the pion-nucleon system, this approach, that was developed 
intensively in the 1990s 0, M, O, inspired a burst of applications m the region of low 
energies [67| . In particular, the time component of the weak axial MECs of the pion-nucleon 
system was constructed [68j within the framework of the heavy baryon chiral perturbation 
theory. Besides the tree approximation, it contains the contribution from the one-loop 
graphs. It was subsequently applied [69| to the calculation of the MECs effect in nuclei that 
has already been discussed above. 

The space component of the weak axial MECs of the pion-nucleon system was constructed 
within the same scheme in Ref. 701 and applied to the weak interaction processes in the 
lightest nuclei in Refs. 

0, 

[7IL l72l. |73|- These calculations are hybrid, since the MECs are 
taken from the EFT, whereas the nuclear wave functions are derived from the potential 
models of the SNPA. Moreover, as it has recently been discussed in [6^ |. the long-range part 
of these MECs does not satisfy the nuclear PCAC constraint. 



Another EFT [66j was applied to the construction of the space component of the weak 
axial MECs in Refs. [74l . u5j where also its influence on the reactions of the low energy 
(anti)neutrinos in deuterons was investigated. 

The validity of the EFTs, based on the nucleonic and pionic degrees of freedom and with 
the heavy meson and A isobar degrees of freedom integrated out, is necessarily restricted 
to the long wave-length limit and internuclear distances r > 0.6 fm [68]. It is natural 
that in parallel with the development and applications of these EFTs, attempts appeared 
76| to construct a class of EFTs without this restriction. The A isobar has already been 



included explicitly within the small scale expansion scheme [65| and it has recently been 
demonstrated in Ref. 77J that under certain conditions one can obtain a consistent power 
counting for model Lagrangians including vector mesons. 

As it is seen from the discussion, a systematic derivation of the weak axial heavy MECs 
operator that would respect chiral symmetry and VDM is lacking. In the absence of an 
EFT including explicitly heavy mesons, we construct here such an operator from our chiral 
Lagrangians in the TAA. Using this operator in calculations of observables at low energies 
and comparing the results with the calculations based on existing EFTs and with the data 
can provide a test of soundness of our approach*. 

The main goal of this study should be seen in the construction of the weak axial nuclear 
exchange currents (WANECs) of the heavy meson range, suitable in the SNPA calculations 
beyond the long wave-length limit, with the nuclear wave functions generated from the 
Schrodinger equation and the related one-boson exchange potentials (OBEPs). For the con- 
struction of the WANECs we make use of the weak axial two-nucleon relativistic amplitudes 
derived in 5(| to which we add the nucleon Born terms. The WANECs are then defined by 
analogy with the electromagnetic MECs 58|, [78j , as the difference between these relativistic 
amplitudes and the first Born iteration of the weak axial one-nucleon current contribution 
to the two-nucleon scattering amplitude satisfying the Lippmann-Schwinger equation. This 

16211 to the construction of the weak axial MECs 



method has already been applied in 



of the pion range. It can be shown in the same manner that the WANECs, defined in this 



* A test of this kind has recently been done in Ref. • 
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way, satisfy the nuclear PCAC equation of the type [79| 



q,Al(2) = [V, + iUmiAl(q 2 )M a (2) , (1.1) 

which follows, as is shortly discussed in appendix |A], from the assumption that the axial 
current consists of one- and two-body terms and it satisfies the PCAC hypothesis for the 
total axial current. Therefore, these currents are suitable for calculations of observables for 
weak processes in the intermediate energy region exactly as the vector exchange currents 
satisfying the nuclear conserved vector current (CVC) equation 

g M y;(2) = [v, v a (i)] , (i.2) 

are applied in the analogous electromagnetic transitions. The most favorable situations 
appear in light nuclei, where the approximation of free nucleons used for construction of the 
transition operators is commonly accepted. In other words, we treat the vector and axial 
currents on equal footing and we consider Eq. (11. ip to be as important for the axial current 
as Eq. 01.21) is important for the vector current. 

The structure of the paper is as follows. In Sect. [Ill we consider the two-nucleon weak 
axial relativistic amplitudes of the p-, u-, and a\ ranges, derived from the Lagrangian 38| 



and we list the PCAC equations which the amplitudes satisfy. In Sect. [Till we first define 
the WANECs as the difference of the relativistic exchange amplitude and of the first Born 
iteration of the nuclear equation of motion. Next we derive the nuclear PCAC that the 
WANECs should satisfy Then we proceed to investigate the structure of the WANECs and 
present the resulting currents in the leading order in 1/M. 

In Sect. HV[ we provide numerical estimates of the cross sections and of the strength of various 
parts of the space component of our WANECs for the weak deuteron disintegration by the 
low energy (anti) neutrinos and compare them with the calculations of Refs. [HI, [75|, [80, 81 



We discuss our results in Sect.|V] Our notations and basic definitions are shortly presented 
in appendix [A], and the structure of the weak axial a meson exchange current is shortly 
discussed in appendix [B] 

II. TWO-NUCLEON WEAK AXIAL MESON EXCHANGE AMPLITUDES OF 
THE p, a x AND u RANGES 

We first write the weak axial amplitude for the ith nucleon (i= 1, 2) 



j^CM) = -<Y<A-v(«)rt,(0 -<Aft M ^(«?)ir = ( 2 -i) 

Here qt =p' ; — p«, the vector-meson propagator is generally designed as 

KM = L v + q -^r) Af( q % A|(g 2 ) = -^—s, (2.2) 

and the pseudovector and pseudoscalar vertices are defined as 

f a 5u (i) = u{ V [) 7 ,7 5 rf tt(ft) , rr = ig fi(pi) irf u( Pl ) = u(p^ 6*r? u{ Vi ) . (2.3) 
The divergence of the amplitude (12.11) is 

q^il, i) = i Um\ AUqf) M a (l, i) . (2.4) 
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Here the one-body pion absorption amplitude M a (l,i) is defined as 

M a (l,i) = V™. (2.5) 
The general structure of the two-nucleon weak axial amplitudes of our model is given in 



Fig.dJ We shall next consider p, a\ and u exchanges. As in the pion exchange case |39|, |40 
we study first the relativistic exchange amplitudes. 

A. Two-nucleon weak axial exchange amplitudes of the p, a\ and uj meson range 

Let us first write down the general form of the nucleon Born amplitude 
Jl,,B = ^b(«i) + ^b(t) = -u{Vx) [6* v) (-q 2 )S F (P) J 5M (l,g)i( G + - a") + J^{l,q) 

u( Pl )A? ) (q 2 )u(p , 2 )di v) (q 2 )u{p 2 ) + (l <- 2) ,(2.6) 



xS F (Q)6? {u) (-q 2 )^a + + a 



where it holds for various meson exchanges 



B 


= P, 




- -i^h - 


B 


= 




= -.y(7, " 


B 




0£ = 





P (Trfiq2s)i = -W P NNO^{q 2 ) , (2.7) 



and 



2M 

- ■^<?r l 5q2s)i = -ig^NNOi V (q 2 ) , T n -»• 1 , (2.8) 

i = -ig ai NNO% , (2.9) 

(2.10) 



a ± = i[7f,7f]±r?. (2.11) 



For the isovector meson exchange 

a + = r 2 a , a' = -i(f ± x r 2 ) a , (2.12) 

whereas for the isoscalar meson exchange 

a + = rf, a" = 0. (2.13) 

As it will become clear soon, the p and a% exchanges should be considered in the chiral 
model 38| simultaneously. 



Other weak axial exchange and pion absorption amplitudes can be derived from the operator 
amplitudes, constructed in Ref. 5(| in conjunction with the Bethe-Salpeter equation, by 



sandwiching them between the Dirac spinors for the free nucleons. The weak axial exchange 
amplitudes of the p meson range are as follows. Besides the nucleon Born amplitude, J^ >p , 
the only potential contact term is J| c „(7r). The mesonic amplitudes belong to the non- 
potential ones and they are aip (ai) and Jg ai (tt) . Considered together with the p 
meson exchange amplitudes, the a\ meson ones contain only the potential amplitudes, that 
are the nucleon Born terms J5L ai and three contact terms Jg Cj ai (it) (i=l,2,3). It can be 
verified that separately, the exchange amplitudes „(2) and ai (2), defined as 

J 5p,p( 2 ) = J 5p,p + J^cpi^) + J5p,a lP i a l) + J^^pi.^) , ( 2 - 14 ) 

6 



and 



■^5/i,ai(2) 



(2.15) 



i=l 



respectively, do not satisfy the standard PCAC constraint. 
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FIG. 1: The general structure of the two-nucleon weak axial amplitudes considered in this paper. 
The weak axial interaction is mediated by the meson B that is either tt or a\ meson. The range 
of the amplitude is given by the meson B2 that is here p, a\ or uj meson; (a),(b) - the nucleon 
Born amplitude ^ b 2 (B); (c) - a contact amplitude c .b 2 (-B); it is connected with the weak 
production amplitude of the B2 meson on the nucleon. Another type of the contact terms, J5 Bi B2 , 
is given in (d), where the weak axial amplitude interacts directly with the mesons B\ and Bi- (e) 
- a mesonic amplitude J5 Bi B2 {B). The associated pion absorption amplitudes correspond to the 
graphs where the weak axial interaction is mediated by the pion, but with the weak interaction 
wavy line removed. There are three types of these amplitudes in our models: Mj| , M® B2 and 
M a BlB2 - 
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However, using a derivation analogous to that of Sect. 3 of Ref. 5(j, one finds § that the sum 
of the amplitudes (12.141) and (12. 15[) satisfies the PCAC equation 



0) + W) =iUmlAl(q 2 ) M a p (2) + M»(2) . (2.16) 



The L) meson exchange amplitudes contain only the potential amplitudes and cw (7r). 
The divergence of the ui meson exchange amplitude Jg w (2), defined as 

yields the PCAC equation 

= i/X A^(g 2 )M«(2). (2.18) 



Using the results obtained in [5Qj, one can derive the two-nucleon weak axial and pion 
absorption amplitudes from the Lagrangian [HlJ in the same manner. We mention only 
that the p and a\ exchange amplitudes satisfy the PCAC separately and that the model 
dependence should be expected at higher energies, because it is of the short range nature. 
In the next section, starting from the obtained two-nucleon weak axial exchange amplitudes, 
we derive the WANECs of the heavy meson range. 



III. WEAK AXIAL NUCLEAR EXCHANGE CURRENTS 

We define the WANEC of the range B as 

i5/Lt,s(2) = J$ti, b(2) — ^5fi, B 5 (3-1) 

where the two-nucleon amplitudes J^ B {2) are derived in the previous section, and t^g 1 is 
the first Born iteration of the one-nucleon current contribution to the two-nucleon scattering 
amplitude, satisfying the Lippmann-Schwinger equation [581 ]. 

P - E(P) + IE 

+ J ? (l,p[,Q)- 1 V B (Q,p^p u p 2 ) + (1 «-> 2), (3.2) 

Q - E(Q) + is 

where P = p[ + q 2 = Pi + q and Q = pi — q 2 = p[ — q, as it is seen from Fig.[2j Further, V B 
is the one-boson exchange nuclear potential and 

3U\p\p) = u{p')J^{l, l)«(p) , (3.3) 

is the weak axial nuclear one-nucleon current in the momentum space. It is the non- 
relativistic reduction of the one-nucleon amplitude J5 (1,1), given in Eq. (12. lft . In the 
nucleon kinematics, q\ = p[ — p\ and q w = E(p[) — E(p\). The amplitude ^^(2) contains 
the nucleon Born terms, graphically presented in Figs.[T](a) and[]Jb). For them, however, 
the four-momentum conservation takes place. Next we derive the continuity equation for 
the WANEC. 



For more details see also Sect. Ill of Ref. [82 
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FIG. 2: The kinematics of the first Born iteration. The nucleon line in the intermediate state is 
on-shell. 

A. The PCAC equation for the WANEC 

In order to derive the nuclear PCAC for the current #(2), we should know the continuity 
equation for the first Born iteration t^J B . This equation is 



5/Lt, 

V B (pi,p^P,p 2 )- m°(l,P,A) 

P - E(P) + xe 



+ m a (l,p[,Q)- r V B (Q,p^]PL,P2) 

Q -E{Q)+ie 



- [V B ,p a 5 (l)\ 

+ (l<-2), (3.4) 

where 

[V B ,p a 5 (l)\ = V B (p{,^;P,p 2 )p a 5 (l,P,pi) - p a 5 (l,p[,Q)V B (Q,p^p u p 2 ). (3.5) 

It follows from Eqs. (12.161) and (12.181) that the divergence of the amplitude B (2) can be 
written as 

q,J^ B (2) = iUmlAl(q 2 )M a B (2) . (3.6) 
Using Eqs. flU} and ([33]), we derive the PCAC for the WANEC, defined in Eq. (jXTft 

Wl B (2) = iUml^ F {f)m B {2) + ([V B ,pi(l)] + (1 <-> 2)) , (3.7) 



where 



m" B (2) = M«(2) -mS FS/ , (3.8) 



and the first Born iteration of the one-nucleon pion absorption amplitude is 



m a B FBI = V B (p[,pi;P,p 2 )- ] , m'Cl.P.fi) 

P - P(P) + ^ 

+ m a (l,#,Q)- r VB(g,^;pi,P2) + (1 <- 2). (3.9) 

Qo - E(Q) + is 
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We shall call the amplitude m B {2) the nuclear two-nucleon pion absorption amplitude of 
the range B. It is seen from Eq. (13. 8p that it is defined in the same way as the WANEC in 
Eq. flEID. 

We now pass to investigate the structure of the WANECs. 



B. The structure of the WANECs 

As it is seen from Eq. (13. ip , the structure of the WANECs differs from the structure of the 
two-nucleon weak axial exchange amplitudes studied in the previous section. In the case of 
the Schrodinger equation, the positive frequency part of the nucleon Born term and the first 
Born iteration differ and this difference provides a contribution to the WANECs. Here we 
shall calculate this difference. For this purpose, we split the nucleon propagator Sf into the 
positive, Sp~\ and negative, Sp~\ frequency parts (see App.|A|). Then the positive frequency 
part J^jj^B °f the nucleon Born term J| B , Eq. (12.61) . can be cast into the form 



Po — & \P) z 



+J 5M (l,p , 1 ,Q) 1 V B (Q,p' 2 ; Pl ,p 2 )l-(a + + a-) + (1 <- 2), (3.10) 
where the quasipotential Vb is defined as 

V B (p' 1 ,p' 2 ;Pl,P2) = 9BNN^{p'l)Ofx){- ( h) U iPl) A (\ V )(<12) U(P2)&( V )(<12)U(P2) , (3.11) 

and Jnfj, is the amplitude Eq. (12. ip without r a /2. For the amplitude t^g 1 , Eq. (13. 2p . one 
can write in a similar way 

ts^B 1 = Vb(p[,P2;P,P2)- — ^ -j^{l,P,p 1 )\{a + - a") 

P - E(P) + 18 * 

+ Ml,p[,Q)~ -V B (Q,p^Pi,P2)l(a + + a-) + (1 - 2),(3.12) 

Q - E(Q) + ie £ 

where Vb and j 5m correspond to Vb and jL, re spectively, but without the isospin dependence. 
Using the definition Eq. (13.11) and Eqs. (I3.10P and (I3.12p . we can calculate the contribution 
to the WANEC arising from the difference of the positive frequency nucleon Born terms and 
the first Born iteration 

A „■<*(+) _ 7 a (+) _ + a,FBI 
a J5/j,B ~~ 5/^, B L 5[i,B 

= ~ Jg.. D^b(P)\^ - a-) + 1 D^ B (Q) I(a+ + a") 

P - E(P) + is 2 Q - E(Q) + ie 2 

+ (1-2), (3.13) 
D,, B {P) = Vfl^p^P,^)^!,^?!)-^^,^;^^)^!,^^), (3-14) 

D^ B {Q) = J5h(1,P'i,Q)Vb(Q,P2,PuP2) - j5^,Pi,Q)VB(Q,P2-,Pl,P2) • (3.15) 



Checking the structure of the currents and potentials, one can see that the contributions to 

:<*(-+ 

5fx, 



Aj'g^g can arise from the difference in the dependence on the energy transfer of the BNN 
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vertices, of the currents, and of the B-meson propagators. We shall call them the vertex, 
external and retardation currents, respectively. They satisfy the equation, 

N a £l = 3l,,B(vert) + j%^ B {ext) + f^ B (ret) . (3.16) 

Besides these corrections, one obtains the whole set of terms by the standard non-relativistic 
reduction of the amplitudes J^^ }B {2) up to the desired order in 1/M. 



C. Results 



In this section, we present the resulting WANECs. 



1. Vertex currents 



As it is seen from Eqs. (12. 7p and (12.81) . the p- and u exchanges can contribute into the vertex 
current j\ B (vert). This contribution can be calculated by expanding in the BNN vertex 

of the pseudopotential V B around Pq = E(P) and Qo = E(Q). Then one obtains 



.B 



V B (p' 1 ,p' 2 ;P,p 2 ) = V B (p[,pi;P,p 2 ) +ig BNN — [Po-E(P)} 
x«(p' 1 )a, 4 n(P)Af ?) (g 2 )^(p 2 )af si (g 2 ) M (p 2 ) 



(3.17) 



and 



.B 



Vb(Q,p' 2 ;pi,P2) = v B (Q,P2,Pi,P2) +Wbnn^[ e (Q) - Qo] 

x u(Q)a J MPi)^Uq2Mp' 2 )0^ st (q2)u(p 2 ) . 
Substituting the expansions Eqs. (13. 17|) and (13 . 18[) into Eqs. (13.13 h — (13.1 51) we obtain 



(3.18) 



jL,B(vert) = ig 2 BNN 



K 



B r 



2M 



u(p' 1 )(Tj4u(P)j5 l ,(l,P,Pi)^(a + - a ) - j 5 »(l,p[,Q) u{Q)a j4 u{pi) 



x -(a + + a ) 
2 V ; 



AfJq 2 )u(j/ 2 )&>' t (< h )u(p 2 ) + (1 - 2) 



(3.19) 



Performing the non-relativistic reduction and preserving only the part ~ (1 + k b ), the space 
part of Eq. (13. 19[) reduces to 



J divert) = ^K B (l + K B )(g A F A { 



(2Mf 
+ ia~ 



a 



q x (a 2 x q 2 ) + i B x x (Pi x (a 2 x q 2 )) 



iPi x (a 2 x q 2 ) - B x x (q x (<j 2 x q 2 )) 



+ 



gp q / . + 



ia 



mi 2M 
q- &x x (q x (<j 2 X q 2 )) 



q- (J\ x (Pi x (a 2 x q 2 )) 



A|(g 2 2 ) 



+ ia 

:i <- 



iq ■ Pi x [a 2 x q 2 ) 
2) , (3.20) 
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where A = p[ + Pi, gp(q 2 )/mi = 2gf 7r Ap(q 2 ), and m l is the lepton mass. This current is 
important for the p meson exchange, since k b = pa 6.6. 

Analogous calculations for the time component jg B (vert) show that in comparison with 
the space component, Eq. (13. 20 p . it is by one order in 1/M suppressed, 



j a 50>B (vert) ^ 0(1/M 4 



(3.21) 



451 ]. only the negative 



According to the definition of the exchange currents in Refs. [25l . 
frequency part of the nucleon Born terms contributes. So the vertex currents are absent 



m 



45 



It means that the currents [45j are suitable for calculations with the nuclear wave 
functions that are solution of the equation of motion, providing the first Born iteration that 
cancels exactly the positive frequency part of the nucleon Born currents, in order to avoid the 
double counting. In the pion exchange current the analogous vertex current of our approach 



62], if the chiral model with the pseudovector 



is again of the nominal order 0{M~ ) [I 
nNN coupling is used. In the chiral model with the pseudoscalar ttNN coupling, the vertex 
current is absent, but the sum of the negative frequ ency part of the nucleon Born term and 



of the PCAC constraint term yields the same result 
included, the pion exchange current [45[ is ~ 0(M~ 



18|. 

-5- 



40, 62 . Without the vertex current 



2. External exchange currents 



These currents arise from the qo dependence of the amplitudes Js M in Eqs. ( 13. 14j) and (13.151) 
One obtains [821 



j% B (ext) 



g A ±A$tf) {-a + [V B 1 \q 2 )(B 1 -(P l + q 2 )) - ■ (A - q 2 ))V^ (£)] 



+a-[V B 1) m(ffi ■ (Pi + &)) + (<?i ■ (A - foMHti)] } + (1 ^ 2) (3.22) 
^ ]A-(g 2 ) - A^g 2 )] { a + [vS\q 2 Wi ■ 0) ~ (?i • fiV$\<h)] 
-a-[VP{q 2 ){a x ■ q) + {a, ■ q)M\q 2 )] } + (1 <- 2) . 



(3.23) 



Here V B is the leading order term in the B-meson exchange potential Vb- In deriving these 
equations, we neglected terms ~ [Af (g 2 )] 2 . 



3. Pair term retardation 



This contribution arises from the pair terms due to the all meson exchanges and it appears 
from the difference between the energy dependence in the propagator of the pseudopotential 
and the potential [5^]. The final result for the retardation currents is 82| 



jl,,B(ret) = ({(1 + u)(P 2 ■ q 2 ) + (1 - v)[{P x + q) ■ q 2 )]}V B 1] (q 2 )] 5 ,(l, P,Pi) 

x \{a + - a') - {(1 + z/)(A • q 2 ) + (1 - v)[{?x - q) • Q) 

x Vg\q2)\{a + + a")) A|(g 2 2 ) + (1 - 2) . (3.24) 
Here v is a parameter of a unitary transformation (cf. [83]). 
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4- Negative frequency part of the nucleon Born terms 



This contribution arises from the nucleon Born currents J§ B , Eq. (12. 6p . by the change 
Sp — > Sp. We consider first the B=p and 00 meson exchanges. 

After the non-relativistic reduction, we obtain for the negative frequency terms of the p- 
and uj ranges 



J5,B 



J 50, B 



^^(1 + k b ) ( a + {g A F A [(at ■ q 2 ) q + q 2 x ((a, + a 2 ) x q 2 ) 



(2M) 

-NPi x (^1 x (ct 2 x 52)) - i(Pi x <?2 
+Pi x ((?! x q 2 



2M (<ti • g 2 ) 9 > 

m ? J 

z(<?i • a 2 x <f 2 )<7 + P\ x (a 2 x <f 2 ) + iq 2 x (<?i x (a 2 x <f 2 )) 



2i M- 



mi 



&i-a 2 x q 2 )q) ) A F (q. 



+ 



9bnn 



a + g A F A [-i(?2 x P 2 ) + Pi x (ai x P 2 



(2M) 3 

- i(JPi x P 2 ) + <f 2 x (cti x P 2 



— a 



{g A F A 



- 2M- 



* ■ P 2 ) q Af (g 2 2 



fa ■ P 2 ) g 
+ (1 432^) 



9bnn 
(2M) 2 



1 + K *) J _ fl H 



g p (q 2 \ 

g A F A i(a 1 ■ &2 x q 2 ) + — g (^i • $2) 



9aFa(vi ■ 92) + z — <?o (*i • 0"2 x q 2 ) 



mi 



- a + g A F A (<7i ■ P 2 ) 



+ a 



0p( 



q (a 1 -P 2 ) A B F {q 



^3.26) 



Comparing our current of Eq. (13. 25|) with the current A±(V) of Ref. [45| we can see a dif- 
ference in the term of the form <7i<f 2 2 . Spurious factor 1/4 appears in [45j due to the use 
of S F ~\p) = («7 • p — M + 74_E'(j» 2 ))/4M 2 , and approximating P = M, when calculating 
the contribution from the 7474 part of the vector exchange. The same is true also for the 
scalar-isoscalar exchange. The correct result is obtained by making use of E = M +p 2 /2M. 
Summing up the currents j£ B (vert), Eq. (13 .20 p . and j^ B \ Eq. (13.251) for B = p, and keeping 
only the terms ~ (1 + kX), we arrive at the nuclear potential term of the p range^ 



V\2 



g A F A {r 2 a [<f x (a 2 x q 2 ) + ia l x (P 1 x (a 2 x q 2 )) 



g P V , (i + " r . 

2 J (2M) 3 

+(fa x r 2 ) a [zP x x (a 2 x <f 2 ) - cti x (q X (a 2 x <f 2 ))] } A£(g 2 



2 7 (2M) 2 
+ (1 <- 2), 



[r 2 (^i ■ g 2 ) - (ri x r 2 ) a (ai • a 2 x q 2 )\ q A p F (q. 2 " 



(3.27) 



™ Here wc follow the nomenclature of Ref. 



62 1. 
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9p 
2 

gp(q 2 



2 (1 + K V ) 

" {i g A F A [t 2 {ff x ■ o 2 x q 2 ) + (t x x f 2 ) a (a 1 ■ q 2 ) 



+ 

mi 

+ (1-2). 



(2M) 2 

<?o [r 2 (<?i ■ 9*2) - (fx x r 2 ) a (ai • a 2 x g 2 )] } A p F (q 2 2 ) 



(3.28) 



The part of our current j£ Qp (pot), Eq. (13.281) . that is ~gA, is in agreement with the current 
J^ a (v—pair), Eq. (38a), derived by Towner 26[. On the other hand, the part of our current 
j£ (pot), Eq. (13.271) . that is ~gA, differs from the p meson pair term, used by Schiavilla et 
ai [H and also in Ref. fl2l |. 

Since Kg ~ —0.12, we keep only the negative frequency Born term contribution into the 
omega meson potential term, 



f) 2 (2P T * ( 9aFa i Pl X X P2) ~ X P2) + (1 + [m ' * 
- z'(Pi x q 2 ) + q 2 x ((ai + ff 2 ) x g 2 ) + iP 1 x (^ x (a 2 x g 2 )) ] } 

(1-2). 



2][ ' JP {l + K S )q(a 1 .q 2 ))A^ 



mi 
9u 



1 



T\ { ■ P-i) + i(! + «s)(*i • o 2 x g 2 )] 



(2M) 2 

+^q (l + Ks )(a 1 -q 2 )} A F (q 2 ) + (1 <- 2) . 



(3.29) 



(3.30) 



In this case, the potential and pair terms coincide. Analogous calculations for the ai meson 
exchange yield the pair term that can be found in Ref. 82J. The currents studied so far are 
derived from the nucleon Born term amplitudes. We now present shortly the currents that 
have the origin in other amplitudes of the section III AI 



5. Potential contact WANECs 



Here we have the currents j^cpi 77 ), 3% 



p., Cia\ 



7T 



i=l,2,3, and jg c£J (7r), obtained by the 



non-relativistic reduction of the pion production amplitudes, entering the weak amplitudes 
^5p, cpl 77 )' J$p, Cicala)' i=l,2,3, and J£ ^(71"), respectively. The explicit results for the currents 
can be found in 82 



6. Non-potential WANECs 

In the considered model, there are two non-potential currents, 7'^, „ ir (ai) and 3$^ aiP (Tt), 
related to the mesonic amplitudes J5 aip (ai) and J^ :aiP (^), respectively [82j. They are of 
the p meson range. In the next section, we shall use the space component of the current 
j% aip (ax) in the numerical estimates of the cross sections. It reads 

+ (a 1 ■ q)(a 2 x q 2 ) + (q 2 ■ a x x o 2 )q{\ A p F (q 2 ) + (1 - 2) . (3.31) 
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We have completed the construction of the heavy meson WANECs of our models. In the 
next section, we apply the leading terms of our currents in calculations of the cross sections 
for the weak deuteron disintegration by the low energy (anti)neutrinos. 



IV. NUMERICAL RESULTS 



We now calculate the cross sections for the reactions of the weak deuteron disintegration by 
low energy (anti)neutrinos, 



v x + d 

V x + d 

v e + d 

V e + d 



v x + n + p, 

K + n + v , 
e~ + p + p, 
e + + n + n , 



(4.1) 
(4.2) 
(4.3) 
(4.4) 



where v x v x refers to any active flavor of the (anti)neutrino. The reactions (14. ip and (14.3P 
are important for studying the solar neutrino oscillations, whereas the reactions (14. 2 p and 
( 14. 4p occur in experiments with reactor antineutrino beams. The precise knowledge of the 
cross sections of the vd reactions is needed [84j for the calculations of the response func- 
tions of the SNO detector (g^, 0, H3, M 
5.21 ± 0.27 (stat.) ± 0.38 (syst.) 
solar models prediction. 

Theoretical studies of the reactions_ (l4. ID - (14.41) . including exchange currents, were accom 



x 



10 6 cm 2 s 



The total active (u x ) 8 Be solar neutrino flux 
1 was found in agreement with the standard 



calculations 80 



plished in theSNPA in Refs. [12 



UHd, and in the EFT's in Refs. 0, H, H- The 



accomplished within the SNPA, generally differ between themselves 

In our opinion, this provides a good motivation 



75 



at the threshold energies by 5%-10% 
to make independent calculations. Before presenting our numerical results, we introduce 
necessary formalism and input data. 



A. Formalism 



Making use of the technique, developed in Refs. [89l . |90| , one obtains the following generic 
equation for the cross sections 



a 



1 



+i 



^M r G 2 w J ^ dx "~" dk' K k' 2 

2 + 23? (. 
hit ^ ' 



+J3J3I < *jf\\Lj\\d> 



£ [j f \<\jf\\Mj\\d> 

^Xj f ,J>0 



j f 3 <\j f \\Mj\\d><Xj f \\Lj\\d> 

2 1 I ^ \ ■ I \mel\ I j ^ |2 



E [l<^/l|rr B IM>r + l<Aj/l|r/||d> 
\>f,J>l 



+ i(jxfh E Z{<\jf\\Tr 9 \\d><\j f \\Tf\\d>* 

Xj f ,J>l 



(4.5) 



Here Gw = Gp for the neutral channel reactions (14. ip . (14.21) . whereas Gw = Gp cos 9q for 
the charged channel reactions (I4.3p . (14.41) . Gp is the Fermi constant and 9c is the Cabibbo 
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angle. Further, the relative momentum of the nucleons in the final state is 



K 



l2M r {k-Ei~A)--\q\ 2 



(4.6) 



where M r is the reduced mass of the nucleons, k(k') is the momentum of the incoming 

(anti)neutrino (outgoing lepton), E[ = \J mf + {k') 2 is the energy of the final lepton. For 
the reactions ( ]4~T]) and ( TO) . 



A = \e d \ 



2.2245MeV, 



whereas for the reaction (14. 3p . 

A = M p - M n + \e d \ = 0.9312MeV r , 

and for the reaction (I4.4p . 

A = M n - M p + \e d \ = 3.5178MeV. 
The 4-momentum transfer is 

q = k — k' , q = k — E[ . 
For the neutral current reactions, Eqs. 04. 1 j) and (14.21) . 

k' max = kx - 4M r + \8M r (2M r -A + k(l- x)) + k 2 (x 2 - V 



(4.7) 
(4.8) 
(4.9) 
(4.10) 

(4.11) 



For the charged channel reactions, Eqs. ( 14. 3 j) and ( 14.41) . the upper limit of the final lepton 
momentum, k' max , is provided by the solution of the equation 



(k - k') 2 + %M r (E[ + A - k) = . 
The lepton form factors entering Eq. 04.5[) are given by the following equations 



J3J3 

Jo jo 

J0J3* 

i{j x f}i 



2[1 - (A;-g)(^-g)] 



4 sin 2 







2[1 - (/34) + 2(M (/?.«?)] 



■[|gl 2 + 2kk'cos 2 - 
Aql cos 2 § 



2[ 1 + (J3 -k)] — > 4 cos 2 



2 ' 



4sg • (k — p) — > — sin - 

1 ^ l?1 2 



9 2 ^ I o . 2 ^ 

g cos - + |g] sin - 



(4.12) 

(4.13) 
(4.14) 

(4.15) 
(4.16) 

(4.17) 



The form factors, presented after the arrows, are valid in the zero mass limit of the outgoing 
lepton. The vector (3 and the unit vectors of the type b are defined as 



P 



k ' 



E, 



1 1 



b 



(4.18) 
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In addition, in Eq. (HTTjl . s = -1(+1) for the reactions (JO) and g^D ( AO) and (Oil ). 
The form factors (I4.13l) - fl4.16l) are the same for all four studied reactions. 
Let us note that for the reaction (14. 3p . the function under integral in Eq. (I4.5P should be 
multiplied by the Fermi function. For the solar neutrinos, a good approximation for this 
function is [75j 

The reduced matrix elements in Eq. (I4.5p read 



< *lA\d]\\d>= — Y, r^ """" ( "' ;rl <(l's)3,\PA\W)l>^, (4.20) 
k Vl Jo r r 



where the nucleon-nucleon partial waves have asymptotics 



<,x^Ui >x sm(K r- l j + 5{), (4.21) 



and the phase shifts and mixing parameters correspond to the Blatt-Biedenharn [91 
vention. 

The multipoles are defined as 



con- 



TZ a9 {q) = ~£-J d^YUq) ■ M , (4-22) 
= / dty ■ j(q), (4.23) 



L lm {q) = -^J- J dn^- X \q) ■ j(q), (4-24) 

Mim = t^- / dn 5 >im(g)jo(g), (4.25) 



where 



^(?) = v^^'^ + v^^ 1 ^' (426) 



Y^iq) = ^^^ m \q)-^^Y^\q). (4.27) 

(4.28) 

The weak hadron neutral current, triggering the reactions (14. II) and (I4.2p . is given by the 
equation 

3nc,» = (1 - 2stn 2 6 w )jl - 2sin 2 6 W ]s^+ 3 3 5fl , (4.29) 

where 6 W is the Weinberg angle [H (sin 2 9 w = 0.23149(15) Q), (j'f ) is the third 
component of the weak vector (axial) current in the isospin space, and js^ is the isoscalar 
vector current. In its turn, the weak hadron charged current is 

fee,, =J:+JL- (4-30) 
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At low energies, the space component of the weak axial hadron current is the most important 
one. 

Our hadron currents consist of the one- and two-nucleon parts. The one-nucleon currents 
are of the form, 

3 a = 2^ i F i p i + iG M^ x Q)\ y , (4-31) 
f 5 = g A F A {a - ^p\Pi° ~ {3-Pi)Pi + (B-q)q-i(P x * <?)]}y • (4-32) 

As to the two-nucleon part, we consider the WANECs only. In addition to the new potential 
exchange currents j£ (pot), j£ w (pair), and the non-potential exchange current j£ aip (cii), 
we include in our calculations the following exchange currents, derived in the chiral invariant 



models [40|, |4J, 162 



1. The 7r potential term, 

x AJ (q 2 2 ) F 2 NN (q 2 ) (a 2 ■ q 2 ) + (1 <- 2) . (4.33) 

2. The p-7r current, 

(7 1 777^ -> 

iCM - -^I'issl'^lW + CH'XI^il)] 

xFpAw^JAJ^JF,™^ 2 ) (? 2 -?2)i(Ti x ? 2 )» + (1 « 2)(.4.34) 

3. The A excitation current of the pion range, 

*°' (A) = SKMA^ FA(?2)( ^ )2[4?2T2 ° + i(<?lXf2)i(fl><?2rl 

xAJffi 2 ) F, 2 K „(* 2 ) • ») + (1 <-> 2) . (4.35) 

4. The A excitation current of the p meson range, 

&p( A ) = - g^f-M) X ^ X ^ T * + ^ lX ^ X ^ X ^ ' ^ X ^ ] 

x A p F (q 2 2 ) F 2 pNN (q 2 ) + (1 «- 2) . (4.36) 



The form factors i* 1 ] , G^-, and Fa, are chosen in accord with Ref. 12j. The deuteron wave 
function and the wave functions of the 1 Sq and 3 Pj states of the final nucleons are obtained 
by solving the Schrodinger equation, making use of the following first and second generation 
realistic potentials: 



(i) The potential OBEPQB |93j, extended to include the Oi meson exchange [94|. The form 
factors, entering the BNN vertices, are of the monopole (dipole) shape for the exchanged 
meson B=7r(p, to, a\). 

(ii) The Nijmegen 93 (Nijm93) and Nijmegen I (Nijml) potentials (95[. In these potentials, 
the exponential strong BNN form factors enter. 
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In order to keep our calculations consistent, we make use of the same BNN form factors 
also in the currents. The value of the itNA coupling, / 2 j VA /47rm 2 r =0.7827 fm -2 , is derived 
from the A isobar width and it is compatible with the pion photo- and electroproduction 
on the nucleon, the value of the constant G\, G*i=2.525, follows from the pion photo- and 
electroproduction (96|, and k^=6.6 j93j . 

In table 1, we present the scattering lengths and the effective ranges, obtained from the 
Nijml, Nijm93, OBEPQG and AV18** [97J potentials, and also the values of these quantities, 
applied in the EFT calculations (75|. For the generation of the final state nucleon-nucleon 
wave functions from the Nijml and Nijm 93 potentials, we adopted the program COCHASE 
(98| . This program solves the Schrodinger equation by employing the fourth-order Runge- 
Kutta method. This can provide the low-energy scattering parameters slightly different 
from those, obtained by the Nijmegen group, which makes use of the modified Numerov 
method [99]. Some refit was necessary, in order to get the required low-energy scattering 
parameters in the neutron-proton and neutron-neutron 1 So states. 

Table 1. Scattering lengths and effective ranges (in fm) for the nucleon-nucleon system in 
the 1 S state, corresponding to the Nijml, Nijm93 Q, OBEPQG 0, AV18 @ potentials 
and as used in the EFT calculations [75], and the experimental values. 





Nijml 


Nijm93 


OBEPQG 


AV18 


EFT 


exp. 


d n p 
J°np 


-23.72 
2.65 


-23.74 
2.68 


-23.74 
2.73 


-23.73 
2.70 


-23.7 
2.70 


-23.740i0.020 1 
2.77 iO.05 1 


Clpp 
r PP 


-7.80 
2.74 


-7.79 
2.71 




-7.82 
2.79 


-7.82 
2.79 


-7.8063±0.0026 2 
2.794±0.014 2 


^nn 


-18.16 
2.80 


-18.11 
2.78 


-18.10 
2.77 


-18.49 
2.84 


-18.5 
2.80 


-18.59±0.40 3 
2.80±0.11 4 


1 Ref. [100J; 2 Ref. 101]; 


3 Ref. [102|; 4 Ref. 1Q3J 



Below, we compare our numerical results for the cross sections with those obtained in the 



pionless EFT [75[ and in the models, developed within the SNPA [12J, l80|, 1811 ] . However, 
in order to make the comparison with the Ref. (75[, one needs to know the constant L\,Ai 
entering the cross sections of the pionless EFT. In the next section, we extract it from the 
numerical values of the cross sections, obtained in various potential models. 



B. Extraction of the low energy constant L\ a 



In Ref. 75J], the effective cross sections for the reactions fl4.1l) - fl4.4p are presented in the form 

veft(E u ) = a{E u ) + L 1>A b(E v ) . (4.37) 



The amplitudes a(E u ) and b(E u ) are tabulated in 75| for each of the reactions (14. lj) — (14.41) . 
from the lowest possible (anti)neutrino energy up to 20 MeV. The const ant L \. a cannot be 



determined from reactions between elementary particles. In our analysis [1041 ] , we extracted 
Li 5 a from our cross sections, calculated in the same approximation, as it was used in Ref. 75j : 
only the 1 5'o wave was taken into account in the nucleon-nucleon final states and the nucleon 



This potential is used in Ref. 80 1 . 
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variables were treated non-relativistically. The knowledge of Li >A allowed us to compare our 
cross sections with <jeft{E v ). The results showed that the EFT cross sections were in better 
agreement with the cross sections of Ref . [80j , even though Nakamura et al. took into account 
also the contribution from the 3 Pj waves and treated the phase space relativistically. We 
defined the extracted value of the constant L^ A as an average value, Li tA , according to the 
equation 

U) , ,.s <? P ot,i — a-i 



Li, 



Li. 



(4.38) 



where (T pa t,i is the cross section, calculated in the potential model and for the i-th 
(anti) neutrino energy. The application of the same equation in the present calculations 
yields the values of L% A , presented in table 2. As it is seen from table 2, the variation of 
L\ t A is about 10 % for the neutral current reactions, whereas it is up to 40 % for the charged 
current reactions. 

Let us note that the values of the constant Li. 4, ob tained from various analyses [l05 . 106], 
are charged with large errors. E.g. , the analysis [106| of the data provided by the experiments 
with the reactor antineutrino beams yielded L\, A — 3.6 ± 5.5 fm 3 . 

Table 2. Values of the constant Li iA (in f mZ ), extracted from the cross sections of the reac- 
tions fj4.lD - fj4.4p . which were calculated with the Nijml, Nijm93, and OBEPQG potentials, 
and from the cross sections NSGK taken from table I of Ref. 18011 . 



reaction 


Nijml 


Nijm93 


OBEPQG 


NSGK 


|4.1| 


4.8 


5.4 


5.0 


5.4 


m> 


5.2 


5.8 


5.4 


5.5 


(I4.3P 


4.4 


5.3 




6.0 


g3D 


4.8 


5.7 


7.2 


5.6 



In the next section, employing the extracted values of the constant L\ }A , we compare our 
cross sections with those of the pionless EFT [75[ , and also with the cross sections of Refs. fl2l . 



80|,|81j 



C. Comparison of the cross sections 

For the comparison presented in table 3 and table 4, we employ the cross sections calculated 
with the nuclear wave functions generated from the Nijml potential. In comparing our 
results with those of Ref. [75[ we make use of the values of the weak interaction constants 
G F = 1.166 x 10~ 5 GeV~ 2 and g A = —1.26, but we adopt the value g A = —1.254 when 



comparing with Refs. [80|, |8l|. In the calculations of the cross sections for the charged 
channel reactions (14.31) and f|4.4p . the value of the Cabibbo angle is taken cos^c = 0.975. 
Table 3. Cross sections and the differences in % between the cross sections for the reactions 
( 14. ip and (14.21) . The first 7 columns is related to the reaction ( 14.11) . In the first column, 
E v [MeV] is the neutrino energy, in the second column, aNijmi (in 10 _42 x cm 2 ) is the 
cross section, calculated with the Nijml nuclear wave functions. Column 3 reports the 
difference between amjmi (Nijml) and the EFT cross section (14.371) <Jeft, calculated with 
the corresponding constant L\ A of table 2. The difference between onsgk-, taken from table 
I of Ref. [80|, and (Jeft, is reported in column 4. Further, A]r 2 ) is the difference between the 
cross sections aNijmi (PNijm93) and unsgk] A 3 is the difference between the cross sections 
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CNijmi and <7yhh, where the cross section oyhh is taken from table I of Ref. [8l| . The second 
part of the table is an analogue for the reaction (14.21) . 







v x + d 




Vx + np 










v x + d 




u' x + np 






E v 


& Nijml 


Nijml NSGK 


Ai 


A 2 


A 3 




°" Nijml 


Nijml NSGK 


Ax 


A 2 


A 3 


3 


0.00335 


1.0 


0.4 


-1.1 


-0.5 


- 


3 


0.00332 


0.3 


0.1 


-1.1 


-0.5 


- 


4 


0.0306 


1.0 


0.2 


-0.8 


-0.2 


12.0 


4 


0.0302 


0.7 


0.2 


-0.8 


-0.1 


9.3 


5 


0.0948 


1.0 


0.2 


-0.9 


-0.2 


5.0 


5 


0.0929 


0.7 


0.1 


-0.8 


-0.1 


0.9 


6 


0.201 


0.8 


0.1 


-1.0 


-0.3 


10.2 


6 


0.196 


0.8 


0.3 


-0.8 


-0.1 


5.8 


7 


0.353 


0.8 


0.1 


-1.1 


-0.3 


8.2 


7 


0.342 


0.5 


0.1 


-0.9 


-0.2 


2.0 


8 


0.552 


0.7 


0.2 


-1.2 


-0.4 


10.1 


8 


0.532 


1.1 


0.8 


-1.0 


-0.2 


3.2 


9 


0.799 


0.8 


0.4 


-1.3 


-0.5 


9.0 


9 


0.766 


0.6 


0.2 


-1.1 


-0.3 


1.0 


10 


1.095 


0.2 


-0.1 


-1.4 


-0.6 


7.8 


10 


1.045 


0.4 


0.2 


-1.2 


-0.4 


-1.5 


11 


1.441 


0.6 


0.5 


-1.6 


-0.8 


9.6 


11 


1.368 


-0.1 


-0.2 


-1.3 


-0.5 


-0.4 


12 


1.836 


-0.2 


-0.3 


-1.8 


-0.9 


8.8 


12 


1.734 


-0.3 


-0.4 


-1.4 


-0.5 


-2.5 


13 


2.282 


-0.1 


0.0 


-1.9 


-1.1 


10.2 


13 


2.144 


-0.3 


-0.2 


-1.5 


-0.6 


-1.7 


14 


2.779 


-0.3 


0.0 


-2.1 


-1.2 


9.9 


14 


2.597 


-0.4 


-0.2 


-1.6 


-0.7 


-3.4 


15 


3.326 


-0.6 


-0.1 


-2.3 


-1.4 


10.8 


15 


3.092 


-0.5 


-0.2 


-1.8 


-0.9 


-3.5 


16 


3.923 


-0.9 


-0.3 


-2.5 


-1.6 


10.5 


16 


3.628 


-0.6 


-0.1 


-1.9 


-1.0 


-4.9 


17 


4.571 


-1.2 


-0.4 


-2.7 


-1.8 


11.2 


17 


4.206 


-0.8 


-0.2 


-2.1 


-1.1 


-5.2 


18 


5.269 


-1.4 


-0.3 


-3.0 


-2.0 


11.1 


18 


4.824 


-1.0 


-0.3 


-2.3 


-1.3 


-6.6 


19 


6.017 


-1.7 


-0.4 


-3.2 


-2.2 


10.6 


19 


5.481 


-1.3 


-0.3 


-2.5 


-1.5 


-6.9 


20 


6.814 


-2.1 


-0.6 


-3.5 


-2.5 


11.7 


20 


6.177 


-1.4 


-0.2 


-2.7 


-1.6 


-8.1 



Comparing this table with table 3 and table 4 of Ref. |l04j we conclude that the effect of the 
3 Pj waves leads our cross sections in better agreement with the EFT cross sections, but the 
effect is weaker in comparison with the one obtained in Ref. (80[. On the other hand, the 
disagreement with the cross sections of Ref. 8lJ is even more pronounced at higher energies. 
Table 3 also shows that our cross sections are in better agreement with the other cross 
sections in the case of the reaction (14.21) . than in the case of the reaction (14.11) . One can 
also conclude from the differences, given in the columns Nijml, NSGK, Ai, and A 2 that the 
cross sections for the reactions (14.11) and (14. 2p are described by both the potential models 
and the pionless EFT with an accuracy better than 3%. 

In table 4, we present the comparison of the cross sections for the reactions in the charged 
channel. 

Table 4. Cross sections and the differences in % between cross sections for the reactions 
L3|) and (14. 4p . For notations, see table 3. 
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v e + d 




P —I— 71 71 

e -|- pp 








v e + d 


— > e H 


—1— 71 71 






°~ Nijml 


Nijml NSGK 


A, 


A„ A.-, 
L-±2 ^3 


F- 


O Nijml 


Nijml NSGK 




A^ Ao 


2 


0.00338 


-5.8 


-0.6 


-7.6 


-6.7 - 






- 


- 






3 


n A/i rr 

0.0455 


-0.8 


-0.3 


-3.0 


-2.0 - 






- 


- 






4 


0.153 


0.1 


-0.6 


-1.9 


-0.9 1.9 






- 


- 






r 

5 


0.340 


1.1 


0.1 


-1.6 


n ^ o n 

-0.6 2.9 


5 


0.0274 


-1.7 


-0.9 


-2.5 


i a m n 

-1.6 9.0 





0.613 


1.5 


0.4 


-1.5 


-0.5 3.1 


6 


0.116 


-0.3 


-0.1 


O "1 

-2.1 


-1.1 8.1 


( 


0.979 


1.6 


0.4 


-1.5 


-0.5 3.1 


t 


n n'7'7 

0.277 


-0.1 


-0.2 


1 O 

-1.8 


-0.7 7.4 


o 

8 


1.440 


-0.4 


-2.4 


-1.6 


-0.6 3.3 


8 


0.514 


0.2 


-0.1 


-1.6 


a r 7 i 

-0.5 7.1 


9 


O AAA 

2.000 


-0.6 


-2.3 


-1.7 


-0.6 3.1 


9 


n o o n 

0.830 


0.1 


-0.2 


-1.6 


-0.5 7.0 


10 


2.662 


-0.2 


-1.7 


-1.9 


n t o o 

-0.7 3.3 


10 


1.226 


0.7 


0.3 


-1.6 


-0.4 6.9 


11 


3.426 


3.4 


3.3 


-2.0 


n n o 1 

-0.9 3.1 


n 


i Tn 1 

1.701 


0.6 


0.2 


-1.6 


-0.5 6.2 


12 


/i on /i 

4.294 


1.0 


0.3 


-2.2 


lion 

-1.1 2.9 


12 


o o r r 

2.255 


0.5 


0.1 


-1.6 


-0.5 6.4 


13 


5.268 


0.8 


0.2 


O A 

-2.4 


1 o o n 

-1.3 2.9 


13 


2.887 


0.3 


0.0 


-1.7 


-0.6 5.9 


1 A 

14 


/■ n /in 

6.348 


0.5 


0.2 


-2.7 


-1.5 2.7 


1 A 

14 


3.596 


0.4 


0.2 


-1.8 


-0.6 5.6 


1 r 

15 


7.535 


0.3 


0.2 


-2.9 


-1.7 2.4 


15 


4.380 


0.1 


0.0 


-2.0 


a 17 r r 

-0.7 5.5 


lb 


s son 


-0.1 


-0.1 


Q 1 
-O. 1 


o n o o 
-z.U 


it) 


K Oil 


0.1 


0.1 


-Z.L 


-u.y o.o 


17 


10.23 


-0.5 


-0.1 


-3.5 


-2.3 1.9 


17 


6.167 


0.0 


0.2 


-2.3 


-1.0 4.4 


18 


11.74 


-0.8 


-0.1 


-3.8 


-2.6 1.1 


18 


7.168 


0.1 


0.4 


-2.4 


-1.1 4.2 


19 


13.36 


-1.0 


-0.0 


-4.1 


-2.9 1.0 


19 


8.237 


-0.2 


0.2 


-2.6 


-1.3 3.9 


20 


15.09 


-1.5 


-0.3 


-4.4 


-3.2 1.0 


20 


9.373 


-0.3 


0.3 


-2.8 


-1.5 3.7 



It can be seen from the left-hand part of table 4 that our cross sections and the cross section 
[80| are smooth functions of the neutrino energy, whereas the EFT cross section exhibits 
abrupt changes in the region 7 < E u < 12 MeV. In our opinion, the reason can be in an 
incorrect treatment of the Coulomb interaction between protons in the EFT calculations. 
It follows from the right-hand part of table 4 that our cross sections and also the cross 
sections by Nakamura et al. 80j are in very good agreement with the EFT cross sections, 



in spite of the fact that the difference up to 3 % between these calculations persists. Large 
difference between our calculations and those by Ying et al. at the threshold energies can be 
understood by a poor description of the neutron-neutron low energy scattering parameters 
by the Paris potential model. 

Next we compare our cross sections with those of Ref. [l2|. For this purpose, we use 
Gp = 1.1803 x 10 -5 GeV -2 and g& = —1.267. The results, displayed in table 5 are for the 
reaction (14. Xp and the OBEPQG potential. The parameters of this potential, required in the 
calculations of the exchange current effects, are g^ NN / ^7r=14.4, ^ JViV -/47r=3.6, g^ NN /4:ir=98, 



A T =8.62 fm 1 , A p =A LlJ =9.38 fm l , and A ai =10.14 fm l . The shape of the table enables one 



to compare our calculations directly with those, presented in table 1 and table 3 of Ref. 12| 
It follows from table 5 that the most important contributions to the cross section are from 
the A excitation currents and from other currents of the pion range. It is also seen that 
the contributions from the ir potential- and the p-ir terms compensate each other to a large 
extent 62|. One can also see that the currents interfere between themselves destructively. 
Besides, the contributions from the heavy meson exchange currents, such as the p potential 
term, the u pair- and a±-p currents, are at the threshold energies numerically insignificant. 
In the last row, we present the cross sections due to the transition 3 S\ — 3 D\ — > 1 Sq only. 
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The comparison of the last two rows shows that the effect of the 3 Pj interaction in the final 
state is « 0.2/0.6/1.2 % at E v = 10/15/20 MeV. 

Table 5. Cumulative contributions to the cross section a u d ( x 10~ 42 cm 2 ) from the weak axial 
exchange currents for various neutrino energies are displayed. The cross section, calculated 
in the impulse approximation (including the relativistic corrections of the order 0(1/M 2 ) is 
presented in the row labelled as IA (5IA). Other contributions correspond to the exchange 
currents as follows: A(tt) -> j 5 3 ^(A); A(p) -> J 5 3 p (A); p(ir) -> Aj 5 a w (pv); p(p) -> j 5 3 p (pot); 

p{oS) — > Js uipo-ir)] p-n — > j^^pn); a^-p — > j5 aip ( a i)- The cross section in the n-th row is 
given by the contribution of all previous currents, the n-th current including. The number 
in the bracket is the ratio of the n-th cross section to the cross section in the row above. The 
first three exchange current contributions are from the long-range exchange currents. The 
cross sections in the last row are obtained with the neutron-proton wave function calculated 
only for the transition 3 Si — 3 D\ — * 1 Sq. 



E v [MeV] 


5 


10 


15 


20 


IA 




0.094 


8(-) 




1.0} 


*8(-) 


3.296 (-) 


6.740 (-) 


+5IA 





0944 


;0.997) 


1 


084 


(0.996) 


3.281 (0.995) 


6.709 (0.995) 


+A(tt) 





0996 


1.055) 


1 


151 


(1.062) 


3.499 (1.067) 


7.175 (1.070) 







0989 


;0.992) 


1 


141 


(0.991) 


3.466 (0.991) 


7.010 (0.990) 


+P-7T 





0997 


1.008) 


1 


152 


(1.010) 


3.502 (1.010) 


7.181 (1.011) 


+A(P) 





0983 


;0.986) 


1 


134 


(0.984) 


3.444 (0.983) 


7.056 (0.983) 


+p(p) 





0986 


1.003) 


1 


137 


(1.003) 


3.455 (1.003) 


7.079 (1.003) 







0986 


1.000) 


1 


138 


(1.001) 


3.456 (1.001) 


7.083 (1.001) 


+O1-/0 





0985 


;0.999) 


1 


136 


(0.999) 


3.452 (0.999) 


7.074 (0.999) 


^(%) 




0.0£ 


184 


1.134 


3.432 


6.993 



However, the currents, adopted in Ref. |12J differ from our currents in several aspects: 



1 . The 7r pair term is constructed from the pseudoscalar ttNN coupling. This makes the 
model incompatible with the chiral invariance. 



2. The p-7i current, used in 



is connected to our current ( 14.341) by the change 
- 2 2ml 



1 + 



m r 



m 2 + q 2 



m 2 + q 2 



(4.39) 



The A excitation currents (I4.35j) and (I4.36P are in [12] supplied with the couplings 
of the constituent quark model, additionally multiplied by a factor 0.8. This factor 
was found to be needed to obtain the experimental value of the Gamow-Teller ma- 
trix element for the triton (3 decay [44| . However, this way of fixing the MECs is 
a model dependent procedure depending on the specific choice of the nuclear model 
which is the AV18 pote ntial in this case. Let us note that the exchange current model 
underestimates |l07| the precise data on the ordinary muon capture in 3 He , 



fi~ + 3 He 



+ 3 H. Besides, the resulting -zrA^A coupling is in a sharp contradiction 



with its value obtained from the A width and also from the pion photo- and elec- 
troproduction on the nucleon. On the other hand, the damping factor M/Ma ~ 0.8 
arises in the vector and axial A excitation currents, if they are constructed from the 
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'gauge symmetric' Lagrangians [108l . 11091 ] . It means that in order to be consistent, 
one should repeat the analysis of the reaction n + p — > d + 7 80] with the damped 



vector A excitation currents and make use of them also in 12 . Let us note that other 



possibility to reconciliate the Gamow-Teller matrix elem ent of the triton f3 decay with 
the data is to vary the cutoff parameters A,,, and A p 
A excitation currents fixed. 



1101 ] , keeping the couplings of the 



4. The structure of our p potential term is different from the p pair term [12j, and the ui 
pair term and the a%-p current are in [l2| absent. 

5. The BNN vertices are of the monopole shape. 

62J, it is not clear which nuclear continuity equation the 



6. As it has been discussed in 
currents [H, 44] satisfy. 



In table 6, we present the results of calculations obtained with our currents modified accord- 
ing to points 2, 3 and 5 above and with the values of the parameters 121 Q^nn 

/4tt=14.81, 

g 2 pNN /4ir=2.0, f^ NA /4nml=0A701 fm~ 2 , A^=4.8 fin" 1 , A p =6.8 fin" 1 . 

In comparing our table 6 with table 3/Model I of Ref. [l2j, we mention our stronger 7r 
potential term and weaker p potential- and p-n terms. Moreover, the tt- and p potential 



terms are of the opposite sign in comparison with the ir- and p pair terms of Ref. [12| . 
The comparison of our p-n row of table 6 with the rows labelled as NPC Bonn and NPC 
AV18 shows that our results are by 1% to 2% smaller, while the same comparison in table 
5 exhibits an opposite effect. 

Table 6. Cumulative contributions to the cross section o v< x (xlO -42 cm 2 ) from the modified 
weak axial exchange currents for various neutrino energies. For the notations, see table 
5. The modification of the currents and the choice of the parameters are described in the 
text. The values of the cross sec tions in the rows labelled as NPC Bonn and NPC AV18 
are displayed for comparison 1 1 if ] . These cross sections correspond to t he calculations [12| . 
restricted to the transition 3 Si- 3 L>i — > 1 Sq, and for the CD-Bonn |l00j ]. and the AV18 [93] 
potentials, respectively. 



E v [MeV] 


5 


10 


15 


20 


IA 


0.0948 (-) 


1.088(-) 


3.296 (-) 




6.741 (-) 


+5IA 


0.0944 (0.997) 


1.084 (0.996) 


3.281 (0.995) 


6 


709 (0.995) 


+A(tt) 


0.0962 (1.019) 


1.107 (1.021) 


3.356 (1.023) 


6 


869 (1.024) 


+pM 


0.0956 (0.994) 


1.099 (0.993) 


3.331 (0.992) 


6 


814 (0.992) 


+ P-7T 


0.0957 (1.001) 


1.100 (1.001) 


3.334 (1.001) 


6 


821 (1.001) 


+A(P) 


0.0955 (0.997) 


1.097 (0.997) 


3.323 (0.997) 


6 


798 (0.997) 


+p(p) 


0.0956 (1.001) 


1.099 (1.002) 


3.329 (1.002) 


6 


810 (1.002) 


+P(w) 


0.0957 (1.001) 


1.099 (1.000) 


3.332 (1.001) 


6 


816 (1.001) 


+Oi-/0 


0.0956 (0.999) 


1.099 (1.000) 


3.330 (1.000) 


6 


813 (1.000) 




0.0956 


1.097 


3.311 


6.732 


NPC Bonn 


0.09552 


1.099 


3.323 




6.765 


NPC AV18 


0.09565 


1.101 


3.332 




6.787 



Comparing table 6 with table 5 one concludes that in this model, the exchange effect is 
much weaker. Now also the p-n current is numerically insignificant. 
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Besides the calculations presented in table 5 and table 6, that are obtained for the reaction 
(14. ip with the OBEPQG wave functions, we calculated also the cumulative cross sections 
for all the reactions (14. ip - (14.41) using the wave functions generated from the Nijml potential 



951 ]. The results are presented in tab les 7 , 8, 9 and 10. In this case, we chose the couplings 



Gpi 9a and cos 9c according to Refs. [112 
Gp = 1.16637(1) x 10~ 5 GeV-\ g A = 



113 



-1.2720(8) 



cos 9, 



0.9730±0.0004±0.0012±0.0002 . 

(4.40) 

Here the Fermi constant Gp is fixed by the muon decay. Making use of this value of Gp, 
the constants g& and cos Qq are extracted from the neutron beta decay. 
Table 7. Cumulative contributions to the cross section <j V( i (xl0~ 42 cm 2 ) for the reaction 
(14. ip . calculated with the nuclear wave functions and couplings and cutoffs, entering the axial 
exchange current operators, that correspond to the Nijml potential, for various neutrino 
energies. For the notations, see table 3. 



c 



E v [MeV] 


5 


10 


15 


20 


IA 


0.0927 (-) 


1.065(-) 


3.227 (-) 


6.599 (-) 


+5IA 


0.0923 (0.996) 


1.061 (0.996) 


3.212 (0.995) 


6.567 (0.995) 


+A(tt) 


0.0977 (1.057) 


1.130 (1.065) 


3.436 (1.070) 


7.047 (1.073) 




0.0969 (0.992) 


1.119 (0.991) 


3.401 (0.990) 


6.972 (0.989) 


+/0-7T 


0.0977 (1.009) 


1.131 (1.010) 


3.439 (1.011) 


7.051 (1.011) 


+A(p) 


0.0962 (0.984) 


1.111 (0.982) 


3.374 (0.981) 


6.912 (0.980) 


+P(P) 


0.0965 (1.003) 


1.114 (1.003) 


3.384 (1.003) 


6.934 (1.003) 


+P(w) 


0.0966 (1.002) 


1.116 (1.002) 


3.392 (1.002) 


6.951 (1.002) 


*CS ) 


0.0966 


1.114 (99.8) 


3.372 (99.4) 


6.868 (98.8) 



Table 8. Cumulative contributions to the cross section a v d (xlO cm ) for the reaction 
( 14. 2p . calculated with the nuclear wave functions and couplings and cutoffs, entering the axial 
exchange current operators, that correspond to the Nijml potential, for various antineutrino 
energies. For the notations, see table 3. 



E p [MeV] 


5 


10 


15 


20 


IA 


0.0908 (-) 


1.016(-) 


2.996 (-) 


5.972 (-) 


+5IA 


0.0905 (0.996) 


1.011 (0.996) 


2.982 (0.995) 


5.942 (0.995) 


+A(tt) 


0.0957 (1.058) 


1.079 (1.067) 


3.198 (1.072) 


6.398 (1.077) 




0.0949 (0.992) 


1.069 (0.990) 


3.165 (0.990) 


6.327 (0.989) 


+/0-7T 


0.0958 (1.009) 


1.080 (1.011) 


3.201 (1.011) 


6.402 (1.012) 


+A(p) 


0.0943 (0.984) 


1.060 (0.982) 


3.138 (0.980) 


6.270 (0.979) 


+P(P) 


0.0945 (1.003) 


1.064 (1.003) 


3.148 (1.003) 


6.291 (1.003) 


+P(w) 


0.0947 (1.002) 


1.066 (1.002) 


3.156 (1.002) 


6.307 (1.003) 


<s ) 


0.0947 


1.064 (99.8) 


3.136 (99.4) 


6.225 (98.7) 



Table 9. Cumulative contributions to the cross section a v d (xlO cm ) for the reaction 
( 14. 3p . calculated with the nuclear wave functions and couplings and cutoffs, entering the axial 
exchange current operators, that correspond to the Nijml potential, for various neutrino 
energies. For the notations, see table 3. 
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E [MeVl 


5 


10 


15 


20 


T A 
1A 


U.oozU {-) 


Z.Doiy-) 


7 onn z' N 

(.oUU (-J 


1 A £0 ( \ 

14. oy {-j 


I XT A 
+01A 


rt oonn /'rt nnv\ 

u.oouy (^u.yy < j 


2.0 1 1 (^u.yyo ) 


< .zoy (^u.yyo j 


14. OZ (u.yyb ) 




U.o4ol (l.UozJ 


von /'i nt".n\ 
z. I Zy (l.Uoy ) 


( . ( oU (l.Ubo J 


1 r ACl I '1 

lo.4y (l.Ubo J 


+p(vr) 


O.o4oo (U.yyo J 


z. ^ oo (o.yyz j 


rbbb (O.yyz ) 


1 r or /n 00 1 \ 

lo.do (o.yyi j 


+/0-7T 


n Q/iQi f 1 mrA 


O 7Q7 M ni 1 
Z. 1 Ol ^l.UHJ 


( . ( OZ ^l.Ull ) 


10.00 ^I.UIZJ 


+A(p) 


0.3438 (0.985) 


2.690 (0.983) 


7.611 (0.982) 


15.23 (0.981) 


+p(p) 


0.3447 (1.003) 


2.698 (1.003) 


7.635 (1.003) 


15.28 (1.003) 


+pM 


0.3452 (1.001) 


2.703 (1.002) 


7.649 (1.002) 


15.31 (1.002) 




0.3451 


2.696 (99.7) 


7.597 (99.3) 


15.12 (99.8) 



Table 10. Cumulative contributions to the cross section o v d (xlO -42 cm 2 ) for the reaction 
(|4.4p . calculated with the nuclear wave functions and couplings and cutoffs, entering the axial 
exchange current operators, that correspond to the Nijml potential, for various antineutrino 
energies. For the notations, see table 3. 



E p [MeV] 


5 


10 


15 


20 


IA 


0.0268 (-) 


1.190(-) 


4.233 (-) 


9.031 (-) 


+5IA 


0.0267 (0.997) 


1.185 (0.996) 


4.213 (0.995) 


8.985 (0.995) 


+A(tt) 


0.0281 (1.051) 


1.259 (1.062) 


4.505 (1.069) 


9.653 (1.074) 




0.0279 (0.993) 


1.249 (0.992) 


4.464 (0.991) 


9.559 (0.990) 


+P-7T 


0.0281 (1.009) 


1.263 (1.011) 


4.518 (1.012) 


9.680 (1.013) 


+A(p) 


0.0277 (0.985) 


1.240 (0.982) 


4.428 (0.980) 


9.475 (0.979) 


+P(P) 


0.0278 (1.002) 


1.244 (1.003) 


4.443 (1.003) 


9.509 (1.004) 


+p(w) 


0.0278 (1.000) 


1.246 (1.002) 


4.452 (1.002) 


9.529 (1.002) 


<S ) 


0.0278 


1.244 (99.8) 


4.426 (99.4) 


9.407 (98.7) 



It is seen from table 7-table 10 that again, the effects of the ir potential- and p-ir terms 

compensate each other considerably. But in this potential model, the p potential- and u 

pair terms interfere additively and the sum of them contribute sensibly. The numbers in the 

brackets show the part (in %) of the 1 Sq cross section from the total cross section. 

In table 11, we compare the 1 S'o phase shifts (in degrees) for the neutron-proton scattering 

obtained from different potentials, used in the calculations of the cross sections, presented 

above. 

Table 11. Comparison of the 1 Sq phas e shi ft (in degrees) for the neutron-proton scattering 
described by the potentials CD-Bonn [iooj . OBEPQG 0, Nijml and Nijm93 @. 



Viab [MeV] 


1 


5 


10 


25 


50 


CD-Bonn 


62.09 


63.67 


60.01 


50.93 


40.45 


OBEPQG 


62.02 


63.56 


59.91 


50.82 


40.13 


Nijml 


62.12 


63.74 


60.10 


51.04 


40.56 


Nij93 


62.05 


63.61 


59.94 


50.85 


40.38 



The phase shifts labelled as CD-Bonn were delivered by R. Machleidt [114j . As it is seen, the 
resulting phase shifts are very close among themselves up to energies of 50 MeV. It means 
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that all four potentials are in the 1 Sq channel about the same for the internucleon distances 
larger than 0.6 fm. 

In Fig.[3j we consider for the reaction (14. II) and for the strongest multipole Tf the ratios of 
the radial density for various exchange currents to the density of the A excitation current 
of the pion range. As it is seen from Eq. (14.51) . the densities depend on the cosine of the 
scattering angle, x = cos 9, and on the energy v' of the outgoing neutrino. In its turn, v' max 
is given by the incident neutrino energy v and x. As it is seen from the figure, the ratios of 
the densities decrease substantially with increasing values of r for the short-range WANECs. 




2 4 6 8 10 

r [fm] 

FIG. 3: The ratios of the radial density of various WANECs to the density of the A excitation 
current of the pion range, A(7r), are presented for the reaction (|4.ip . and for the Tf 1 multipole. 
The selected kinematics is v = 5 MeV, 8 = 45°, v' = v' max /2. With this choice, the momentum 
transfer is |gj = 4.12 MeV. The curves correspond to the following WANECs: 1- A(p), 2- p(7r), 
3- p(p), 4- p(u;), 5- p — ir, 6- a\ — p. The sign + (-) at a branch of a curve means that the ratio 
of the densities is positive (negative). The signs of the branches before and/or after the kinks are 
opposite. 
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V. DISCUSSION OF OBTAINED RESULTS 



In this paper, we construct the weak axial nuclear heavy meson exchange currents of the p-, 
oj- and ai meson range in the TAA. These currents are suitable for calculations employing 
the nuclear wave functions obtained by solving the Schrodinger equation with the OBEPs. 
Including the ir exchange current, our model currents provide the realistic space compo- 
nent of the axial exchange currents. Adopting for the exchange currents the needed input 
(couplings, cutoffs) from the OBEPs, we can made a reliable estimate of the effect due to 
these exchanges. Since our Lagrangians belong generically to the non-linear realization of 
the chiral group, the a meson is explicitly lacking. On the other hand, the OBEPs contain 
the a meson exchange, needed to explain the medium range attraction. Then it follows 
that the exchange effect, arising from this meson, cannot be calculated consistently in our 
approach. However, the structure of the space component of the WANECs of the o range 
does not allow it to play any significant role at low energies. This is discussed in more detail 
in appendix [Bj 

In accord with Ref. [581, we define the nuclear exchange currents in Eq. fl3.ll) as the dif- 
ference of the relativistic Feynman amplitudes of Sect. [Ill and of the first Born iteration of 
the nuclear equation of motion (J372]). The currents, derived in this manner, are given by 
a finite set of terms and they satisfy the nuclear PCAC equation (13.71) . Let us note that 
our approach follows methods standardly used in the few-body nuclear physics at low and 
intermediate energies: the transition operator that is constructed for free nucleons is sand- 
wiched between the nuclear wave functions obtained by solving the Schrodinger equation 
with realistic potentials. As is well known from the textbooks, if the transition operator 
is the sum of the one- and two-nucleon vector currents, it should satisfy the nuclear form 
of the CVC hypothesis (11.21) . Then accepting chiral symmetry as one of the corner-stones 
on which the nuclear physics calculations rely, one should construct the axial currents sat- 
isfying the nuclear PCAC constraint, as it is done in this paper. If the particles cannot be 
considered as free (nuclear matter) or non-relativistic then one should use other methods to 
describe the nuclear states, and to construct the currents in conjunction with the relevant 
equation of motion for nuclear states, as it w as done, e.g., in Ref. 5(| for the axial currents 
and the Bethe-Salpeter equation, or in Refs. [m, [H3 for the electromagnetic currents and 
the Blankenbecler-Sugar-Logunov-Tavkhelidze quasipotential equation. 
After performing the standard non-relativistic reduction, we obtain the currents in the 
leading order in 1/M. The non-zero difference arises due to the difference between the 
positive frequency part of the nucleon Born terms and of the first Born iteration, which 
gives rise to the vertex, external and retardation terms. The sum of the vertex terms and 
of the negative frequency Born terms (pair terms) results in the nuclear potential currents. 
The vertex currents contribute essentially in the case of the 7T- and p meson exchanges. Note 
that the space parts of our n- and p potential currents (I3.27P differ from the the space parts 



of the 7i- and p pair terms, used in Refs. [12j, |44j, |45 



Besides the contribution from the nucleon Born terms, the potential contact and non- 
potential currents appear after the non-relativistic reduction. These currents have never 
been constructed before. An infinite set of such terms was supposed to exist in Ref. 451 ]. 
where it was argued that these terms should be numerically insignificant because of the 
short-range suppression present in nuclear wave functions. We have shown here that the 
method based on the chiral invariance provides a finite set of terms, of which contributions 
can be estimated quantitatively. 
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The application of our currents to the reaction of weak deuteron disintegration by 
(anti)neutrinos at low energies shows that our currents differ from the set that has recently 
been used in Ref . [12| , in that (i) The tt- and p potential terms contribute with the opposite 
sign, in comparison with the 7r- and p pair terms [121 ]. respectively, (ii) The effect of the 
7r potential term cancels strongly the effect, arising from the p-n current. Besides we have 
shown that at the threshold energies, the contributions to the cross sections from the heavy 
meson exchanges interfere destructively in conjunction with the OBEPQG potential, and 
the heavy meson effect is strongly suppressed in comparison with the effect arising from the 
currents of the pion range, as supposed in [451 ] . However, in conjunction with the potential 
Nijml, the p potential- and to pair terms interfere additively, and this sum is non-negligible. 
This result is important in particular in the case, if the A isobar currents are adopted with 
the suppressed strength [l2j. Such analysis can be performed for other potentials and/or 
energies, too, thus having bookkeeping of the exchange current effects under control. But 
at higher energies, other parts of our currents, not considered in these calculations, should 
be also estimated. 



We compared our cross sections with those obtained in the EFT [74], |75J in terms of the 
constant L\^a- Our values of the constant L\ t A are somehow smaller than its value obtained 
in the analysis of Refs. [751. 80) , but in a good agreement with the dimensional analysis and 



with the data analysis [l05l ~ 106 ] . The comparison of the cross sections for the reactions 



in the neutral current channel and for the reaction (14.41) shows that they are described by 
both the potentials models and the pionless EFT with an accuracy better than 3 %. For 
the reaction (14.31) . the achieved accuracy is ~ 3.4 %. In this case, abrupt changes in the 
EFT cross section in the interval 7 < E u < 12 MeV appeared, whereas the cross sections 
obtained from the potential models are smooth. 

Our exchange charge and current densities are obtained in a general reference frame in terms 
of individual nucleon coordinates. This is a good approximation in the considered energy 
region, since the estimated effect of the boost currents is negligible. For higher energies, 
the calculations of the observables by sandwiching the operators between the intrinsic wave 
functions require t o go one step further: by analogy with the case of the electromagnetic 



exchange currents [117] , the center-of-mass frame dependence should be separated. This 



will be done elsewhere. 
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APPENDIX A: DEFINITIONS AND NOTATIONS 



We use the Pauli metrics. The free nucleon spinors are normalized according to 

u + {p)u{p) = v + {p)v{p) = 1 . 

The spinors are 



(Al) 



u(p) 



E(p) + M 



1 



v(p) 



E(p) + M 



<T-p 



E(f)+M 
1 



(A2) 



E{p)+M 

The nucleon propagator is split into the positive- and negative frequency parts as follows 



where 



S { F +) (P) 



Sf( P ) 



1 



1 



s { ;\ P ) + sp(p) 



+ M 



H/ 



-u(p)u(p), sy(p) 



Po - E(p) p + E(p) 

For a weak semileptonic reaction in the two-nucleon system NN, 



v (~ p) v (— p) ■ 



NN(Pi) + k{ Pi ) — NN(P f ) + l f {p f ), 
we write the field-theoretical S-matrix element in the form 

S = l {2^5^\P f + Vf -P i - Vi )l^)W a M), 
where the matrix element of the lepton weak current is 

Z M (0) = (lf,P f \U0)\h,Pi) , 
the weak hadron current consists of the weak vector and weak axial vector parts, 

W;{q) = J;{q) + J a 5fl (q) , 



(A3) 



(A4) 



(A5) 



(A6) 



(A7) 



(A8) 



and the momentum transfer q = pi — pj = Pf — Pi. In what follows, we deal only with the 
weak axial hadron current J^{q)- 

We define the potentials and quasipotentials in a similar way, as it is done in Appendix A 



of Ref. 58 



We discuss next how the operator of a weak axial nuclear current j^, used in conjunction 
with the equation describing our NN system, can be related to the field-theoretical current 
Jg M (g). Let the time evolution of the NN system be described by a Hamiltonian H = T + V, 
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where T is the kinetic energy and V is the nuclear potential, which we take as the sum of 
the OBEPs, 

V = ]T V B . (A9) 

B=TT,p,uj,ai... 

If the nuclear physics calculations, based on the current and on the eigenfunctions of 
the Hamiltonian H should reflect the PCAC, then the current should satisfy the continuity 
equation 

q-fM = [H, p\{q)] + if«mlAW)™ a (q), (A10) 

where m a is the associated pion absorption amplitude. Supposing that the current consists of 
the one- and two-nucleon components, this equation splits into the following set of equations 

qi-J a 5 (l,q}) = [T l ,p%(l,q 1 )] + if 7T mlAUq 2 )m a (l,q i ), 2 = 1,2, (All) 
q-f 5 (2,q) = [T 1 + T 2 ,p a 5 (2,q)] + ([V , p a 5 (l, <?)} + (1 <- 2)) 

+ if w mlAUq 2 )m a &q). (A12) 

In Eq. (1A12|) . we neglected p%(2, q) in the second commutator on the right hand side. Tak- 
ing into account that q = q w + q 2 o, we find that Eqs. (13.71) and flA12j) are in the full 
correspondence in the space of the nuclear states, that are described by the eigenfunctions 
of the Hamiltonian H. So we can consider the current j\ s(2), defined in Eq. (13. ip . as the 
WANEC of the range B. 



APPENDIX B: THE WANECS OF THE a MESON RANGE 

Realistic OBEPs contain standardly the a meson exchange, describing the attraction at 
medium distances. On the other hand, our WANECs do not contain the component due 
to this exchange so far. This is due to the fact that our Lagrangians reflect the non-linear 
chiral symmetry and it is not clear how to include the a meson into the scheme consistently. 
Here we construct the axial exchange current of the a range starting from a Lagrangian 

AC = 9(T NN^ + i^g a N lb {f.jf)N(j) + ig wNN N l5 (T-ir)N + ig A g p N lul5 (r ■ a v )N , (Bl) 

where besides the pseudoscalar tcNN coupling standardly accepted aNN coupling is present. 
The relativistic amplitudes derived from this Lagrangian by analogy with Sect. [IT] for other 
exchanges of our model are the nucleon Born amplitude J| a and the only potential contact 
term Jg co .(7r). It can be verified that these amplitudes satisfy the PCAC equation 



JL. + JLe. =iUmlAl{q 2 )[M a a + AC]. (B2) 



Using the methods developed in Sect. lIIIt we obtain the pair term of the a range 
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Impair) = ig A F A j0^ [q 2 x (Pi + q) + iP x x {a, x q)} A^(g 2 2 ) + (1 <- 2XB3) 

J^Apair) = 9aF A j^ (^ • (P 1 + q)) t« A° F (q 2 2 ) + (1^2). (B4) 

For the transition 3 Si — 3 Di —* 1 Sq, only the second term at the right hand side of the space 
component of the current j 5 a a (pair) contributes. However, being proportional to q, it is of 

little importance at the threshold. In contrast to the current A±(S) presented in Eq. (2.5a) 



of the Ref. 45J, our current (1B3P does not contain the spurious term analogous to — a\k 2 /4. 
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